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CHAPTER  11 
THE  DATA  SETS 


The  methods  developed  in  the  previous  chapters  were  applied  to  twelve  differ¬ 
ent  Bets  of  data.  These  will  now  be  described. 

11.1  Primary  Mental  Abilities 

Thurstone  and  Thur stone  (1941)  administered  6o  tests  to  710  eighth  grade 
students.  The  intercorrelations  analyzed  were  taken  from  nine  of  these  tests. 

The  first  three  of  these  were  verbal  tests.,  the  next  three  spatial,  and  the  last 
three  numerical,  as  follows: 

1.  Sentences 

C.  Vocabulary 

3*  Completion 

4.  Flags 

5.  Figures 

6.  Cards 

7.  Addition 

8.  Multiplication 

9*  Three  higher 

11.2  Twenty -Four  Psychological  Tests 

This  set  of  correlations  comes  from  a  battery  of  twenty-four  psychological 
tests  given  to  145  seventh  and  eighth  grade  school  children  in  a  subtirb  of  Chicago. 
The  initial  data  were  gathered  by  Holzinger  and  Swineford  (1939)*  The  data  have 
subsequently  been  analyzed  by  a  number  of  investigators  including  Holzinger  and 
Harmon  (1941),  Kaiser  (1958),  Neuhaus  and  Wrigley  (1954),  Harmon  (1967),  and 
others,  so  that  the  characteristics  of  the  data  have  come  to  be  well  known.  The 
tests  are  identified  as  follows: 
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1. 

Visual  perception 

13- 

Straight -curved  capitals 

2. 

Cubes 

14. 

Word  recognition 

3- 

Paper  form  board 

15. 

Number  recognition 

k. 

Flags 

16. 

Figure  recognition 

5- 

General  information 

17- 

Object -number 

6. 

Paragraph  comprehension 

18. 

Number -figure 

7- 

Sentence  completion 

19. 

Figure -word 

8. 

Word  classification 

20. 

Deduction 

9- 

Word  meaning 

21. 

Numerical  puzzles 

10. 

Addition 

22. 

Problem  reasoning 

11. 

Code 

23- 

Series  completion 

12. 

Counting 

24. 

Arithmetic  problems 

11.3  Thirty -Three  Variable  Speed  Study 

These  data  are  from  a  study  by  Lord  (195&)  designed  to  investigate  the  speed 
factor.  Tests  were  administered  to  64 9  students  in  the  entering  class  at  the  United 
States  Naval  Academy  at  Annapolis.  The  tests  were  designed  to  measure  verbal,  ~ 
spatial,  and  arithmetic  reasoning  ability.  In  each  area,  seven  tests  were  admin¬ 
istered.  One  was  the  regular  admissions  examination  denoted  by  (a).  The  remaining 
six  were  short  experimental  tests  parallel  in  content  but  different  in  degree  of 
speededness.  Two  designated  (L)  involved  virtually  no  speed,  one  was  moderately 
speeded  (M),  and  the  remaining  three  (s)  were  highly  speeded.  Six  reference  facto- 
tests  designated  by  (R)  were  also  included.  In  addition,  grades  in  six  areas 
designated  (G)  were  included  as  variables.  The  33  variables  are  as  follows: 
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1. 

Word  fluency  (r) 

18. 

Arithmetic  reasoning  (L) 

2. 

Verbal  (A) 

19- 

Arithmetic  reasoning  (M) 

3* 

Vocabulary  (L) 

20. 

Arithmetic  reasoning  (S) 

4. 

Vocabulary  (L) 

21. 

Arithmetic  reasoning  (S) 

5- 

Vocabulary  (M) 

22. 

Arithmetic  reasoning  (S) 

6. 

Vocabulary  (s) 

23* 

Number  speed  (R) 

7- 

Vocabulary  (S) 

24. 

Number  speed  (R) 

8. 

Vocabulary  (S) 

25. 

Cancellation  (R) 

9- 

Spatial  relations  (A) 

26. 

Picture  discrimination  (R) 

10. 

Intersections  (l) 

27. 

Number  checking  (R) 

11. 

Intersections  (L) 

28. 

English  (G) 

12. 

Intersections  (M) 

29. 

Foreign  language  (G) 

13. 

Intersections  (S) 

30. 

Engineering  drawing  and 

descriptive  geometry  (G) 

14. 

Intersections  (S) 

31. 

Chemistry  (G) 

15. 

Intersections  (S) 

32. 

Mathematics  (G) 

16. 

Mathematics  (A) 

33- 

Conduct 

17. 

Arithmetic  reasoning  (L) 

11.4  Thurstone  Twenty -Variable  Box  Problem 

These  data  are  from  the  classical  study  by  Thurstone  (1947)  designed  to 
illustrate  the  principle  of  simple  structure.  Measurements  of  a  random  collection 
of  thirty  boxes  were  made.  The  three  dimensions  X,  Y,  and  Z  were  recorded  for 
each  box.  A  list  of  2 6  arbitrary  score  functions  was  then  prepared.  Twenty  of 
these  functions  were  included  as  variables  in  our  analysis.  These  are  as  follows: 


1.  X 


11.  A 

2.  Y  12.  TLd 


3. 

z 

13. 

2X  +  2Y 

4. 

XY 

14. 

2X  +  2Z 

5. 

xz 

15. 

2Y  +  2Z 

6. 

YZ 

16. 

/X2  +  Y2 

7. 

y?Y 

17. 

/x2  +  z2 

8. 

XY2 

18. 

/y2  +  z2 

9. 

x2z 

19. 

XYZ 

10. 

xz2 

20. 

2  £  £ 
/X  +  Y  +  Z 

11.5  Eight -Variable  Body  Typ*  Measures 

These  data  are  from  a  study  of  eight  physical  variables  by  Mullen  (1939) • 
The  data  have  been  used  for  illustrative  purposes  by  Harmon  (1967)  and  by  Kaiser 
and  Caffrey  (1965).  The  variables  are  as  follows: 


1. 

Height 

5*  Weight 

2. 

Arm  span 

6.  Bitrochanteric  diameter 

3- 

Length  of  forearm 

7.  Chest  girth 

4. 

Length  of  lower  leg 

8.  Chest  width 

11.6  Twelve -Variable  Anthropometric  Measures 

These  data  are  from  a  factor  analysis  by  Hammond  (1942)  involving  twelve  body 


measurements  on  adult  men.  Hammond  attempted  to  interpret  the  resulting  factor 
matrix  without  rotation  of  axes.  Later  Thurstone  (1946)  reanalyzed  the  data 
rotating  to  simple  structure.  The  variables  are  as  follows: 
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1. 

Stature 

7. 

Chest  depth 

2. 

Sitting  height 

8. 

Head  length 

3. 

Shoulder  breadth 

9- 

Head  breadth 

k . 

Hip  breadth 

10. 

Head  height 

5- 

Span 

11. 

Hand  length 

6. 

Chest  breadth 

12. 

Hand  breadth 

11. T  Fifteen  Variables  from  Hemmerle 

These  data  are  from  a  study  by  Hemmerle  (1965)  designed  to  illustrate  a  method 
for  obtaining  maximum  likelihood  estimates  of  factor  loadings  and  communalities 
using  an  iterative  computer  procedure.  Later  the  data  were  reanalyzed  by  methods 
developed  by  Joreskog  (1967)  and  by  Horst  (1968b).  This  data  set  was  included 
because  of  the  divergent  results  obtained  by  the  several  investigators.  Hemmerle 
does  not  indicate  the  source  of  the  data,  the  number  of  cases,  nor  the  nature  of 
the  variables. 

11.8  Seventeen -Variable  Data  fran  Be cht old— Sample  1 

These  data  are  from  a  study  by  Bechtold  (1961)  designed  to  investigate  the 
factor  analysis  stability  hypothesis.  The  data  are  a  portion  of  tho~~  originally 
collected  by  Thur atone  and  Thurstone  ( 19^1 ) •  The  study  included  seventeen 
variables  from  a  sample  of  212  cases.  The  first  tvo  variables  were  designed  to 
measure  memory  (M),  the  next  three  verbal  ability  (V),  and  successive  sets  of 
three  measure  word  fluency  (W),  spatial  ability  (s),  number  ability  (N),  and 
reasoning  ability  (R).  The  seventeen  variables  were  given  designations  as 


follows: 
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1. 

First  names  (M) 

10. 

Figures  (S) 

2. 

Word -number  (M) 

11. 

Cards  (S) 

3- 

Sentences  (V) 

12. 

Addition  (N) 

4. 

Vocabulary  (v) 

13. 

Multiplication  (N) 

5- 

Completion  (V) 

14. 

Three  higher  (ll) 

6. 

First  letters  (W) 

15. 

Letter  series  (R) 

7. 

Four -letter  words  (W) 

16. 

Pedigrees  (R) 

8. 

Suffixes  (W) 

17. 

Letter  groupings;  (R) 

9- 

Flags  (S) 

11.9  Seventeen -Variable  Data  from  Bechtold- -Sample  2 

These  data  are  from  the  same  study  by  Bechtold  (1961)  as  those  in  Section 
11.8.  The  variables  are  the  same  as  in  that  data  set  but  the  cases  are  a  separate 
sample  of  213  cases.  The  two  samples  of  cases  were  formed  by  assigning  each  of 
425  cases  alternately  to  one  or  the  other  of  two  groups  after  the  cases  were 
thoroughly  randomized. 

11.10  Nine -Variable  Synthetic  Data 

The  correlation  matrix  for  this  data  set  was  derived  from  a  configuration  of 
points  constructed  so  as  to  provide  a  severe  test  for  the  simple  structure  trans¬ 
formation  procedure  described  in  Chapter  9.  A  right  spherical  triangle  was  con¬ 
structed  on  the  surface  of  a  sphere  of  unit  radius.  A  point  was  located  on  each 
side  of  the  triangle  midway  between  the  two  vertices,  or  45  degrees  from  each  of 
the  two  vertices.  Two  more  points  were  located  on  each  side  of  the  spherical 
triangle,  one  each  midway  between  a  vertex  and  the  mid  point,  or  22.5  degrees  frm 
a  vertex  and  the  mid  point.  Thus  the  three  points  on  the  90  degree  arc  of  the 
great  circle  constituting  a  aide  of  the  triangle  divided  the  side  into  four  eoia; 
arcs  of  22.5  degrees  each.  The  cosines  of  the  angular  distances  between  all  pairt 
of  the  nine  points  were  calculated  to  obtain  a  correlation  matrix.  The  ccsinei; 
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of  the  angular  distances  of  each  of  the  nine  points  with  each  of  the  three  vertices 
of  the  right  spherical  triangle  were  also  calculated.  These  values  are  the  simple 
structure  factor  loadings  of  the  variables  (points).  An  adequate  method  of  analysis 
of  the  correlation  matrix  including  simple  structure  transformation  should  recover 
the  synthetically  constructed  simple  structure  factor  leadings. 

11.11  Reading  Comprehension  Factors 

These  data  are  from  a  study  by  Davis  (1944)  designed  to  investigate  the  prim* 
ary  factors  of  reading  comprehension.  Tests  were  designed  to  measure  nine  different 
reading  skills.  The  correlations  are  based  on  scores  of  421  college  freshman.  711' - 
tests  were  as  follows: 

1.  Knowledge  of  word  meaning 

2.  Ability  to  select  the  appropriate  meaning  for  a  word  or  phrase  in 
the  light  of  its  particular  contextual  setting 

3«  Ability  to  follow  the  organization  of  a  passage  and  to  identify'  ante¬ 
cedents  and  references  to  it 

4.  Ability  to  select  the  main  thought  in  a  passage 

5»  Ability  to  answer  questions  that  are  specifically  answered  in  a 
passage 

6.  Ability  to  answer  questions  that  are  answered  in  a  passage  but 
not  in  words  in  which  the  question  is  asked 

7.  Ability  to  draw  inferences  fro®  a  passage  about  its  contents 

6.  Ability  to  recognize  the  literary  devices  used  in  a  passage  and  tr 
determine  its  tone  and  mode 

9*  Ability  to  determine  a  writer’s  purpor  intent,  and  point  of  view, 
i.e. ,  to  draw  inferences  about  a  writer. 
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11.12  The  Heyvood  Case 

These  data  are  from  a  five -variable  synthetic  example  from  Thomson  (1950). 

The  correlation  matrix  was  constructed  so  that  every  tetrad  difference  is  exactly 
zero,  but  the  g  factor  saturation  for  one  of  the  tests  is  greater  than  unity.  How¬ 
ever,  the  matrix  is  positive  definite.  This  example  was  included  to  test  the 
behavior  of  various  scaling  and  loss  function  parameters  described  in  Chapter  8 
for  the  Heyvood  case.  The  loadings  of  the  variables  for  the  g  factor  were  chosen 
as  follows: 


1. 

1.05 

2. 

•  9 

3. 

.8 

4. 

.7 

5. 

.6 

V 


CHAPTER  12 

EXPERIMENTAL  RESULTS 

In  this  chapter  ve  shall  merely  present  the  numerical  results  of  the  analyses 
for  the  twelve  data  sets.  In  Chapter  13  we  shall,  discuss  some  of  the  more  interest¬ 
ing  of  these  results.  At  the  end  of  this  chapter  the  results  are  presented  succes¬ 
sively  for  each  of  the  twelve  data  sets.  For  each  data  set  six  separate  sets  rf 
analyses  are  presented.  For  the  first  group  of  three  of  these  analyses  the  less 
function  parameter  P^  =  1  was  used  and  for  the  second  group  of  three  this  parameter 
was  *  0.  Within  each  group,  the  fi.  st  se  is  for  the  scaling  parameter  p  *  0, 
the  second  for  p  =  .5,  and  the  third  for  t>  1.0.  The  format  for  all  sets  of  data 
is  identical.  It  consists  of  a  first  line,  a  second  line  or  sequence  of  lines,  a 
third  block  of  lines,  and  a  final  line.  These  we  shall  now  interpret. 

12.1  The  First  Line 

The  first  line  has  three  successive  groups  of  numbers.  The  first  group  con¬ 
sists  oy  six  integers.  The  second  consists  of  three  figures.  The  third  group 
has  figures  equal  in  number  to  the  number  of  roots  m  of  the  correlation  matrix 
greater  than  unity. 

The  six  integers  in  the  first  group  are  as  follows: 

(1)  The  first  integer  is  simply  the  arbitrary  serial  order  of  the  data  set. 

(2)  The  second  integer  is  the  order  of  the  correlation  matrix  or  the  number 
of  variables  n. 

(3)  The  third  integer  is  the  number  of  factors  solved  for.  This  is  the 
number  of  roots  of  the  correlation  matrix  greater  than  unity.  It  is  the  same  as 
the  number  of  figures  in  the  last  group  of  the  first  line. 

(4)  The  fourth  integer  is  a  code  for  the  loss  function  parameter  P^.  For 

P  *  1  the  integer  is  1  and  for  P„  *  0  the  integer  is  2.  Thus  th.  integer  1  means 
W  W 

that  the  loss  function  includes  only  the  residual  covariance  elements,  while  the 
integer  2  means  that  it  includes  both  the  unit  weighted  residual  variance  and 

'  . .  **  '  ■*  - - -  •  "  '  • . 


covariance  elements.  It  is  possible,  of  course,  to  have  take  any  value  between 
unity  and  zero  but  only  the  two  extremes  were  used  far  all  twelve  data  sets. 

(5)  The  fifth  integer  is  a  code  for  the  scaling  function  parameter  p.  For 
p  =  0  the  integer  is  1,  for  p  =  .5  it  is  2,  and  for  p  =  1.0  it  is  3*  Thus  the 
integer  1  means  that  the  scaling  function  is  the  square  root  reciprocal  of  the 
residual  variance,  the  number  2  means  that  it  is  the  square  root  residual  of  the 
total  variance,  and  the  number  3  means  that  it  is  the  square  root  residual  of  the 
estimated  or  common  variance.  It  is  possible,  of  course,  to  let  p  take  any  value 
between  zero  and  unity  but  only  the  three  values  indicated  above  were  used  for  the 
twelve  data  sets. 

(6)  The  sixth  and  final  Integer  in  this  group  is  a  code  to  indicate  the  rcv 
scaling  treatment  of  the  factor  loading  matrix  prior  to  the  simple  structure  pro¬ 
cedures  of  Chapter  9-  The  integer  1  indicates  that  the  factor  loading  matrix  was 
normalized  by  rows  prior  to  simple  structure  transformation;  if  it  was  not,  the 
integer  2  is  used  to  so  indicate.  The  computer  program  provides  for  both  options 
but  in  this  study  only  the  normalizing  option  was  used  for  all  sets  of  data.  Hence 
for  each  of  the  six  analyses  for  all  twelve  data  sets,  the  last  integer  in  the 
first  group  of  six  in  the  first  row  is  always  1. 

The  three  figures  in  the  second  group  are  as  follows: 

(1)  The  first  figure  in  this  group  is  the  ratio  of  the  sum  of  squares  of  the 
first  m  roots  of  the  matrix  for  specified  scaling  and  loss  function  parameters  to 
the  sum  of  squares  of  all  the  roots  of  this  matrix.  This  is  the  criterion 
developed  in  Chapter  8  which  it  is  desired  to  maximize.  The  maximum  value  it  can 
attain  is  unity. 

(2)  The  second  figure  in  this  group  is  the  number  of  iterations  required  to 
reach  the  tolerance  limits  for  the  equality  of  two  successive  iterations  for  0  or 
the  iteration  limit,  whichever  is  reached  first. 
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.  (3)  This  figure  is  the  time  in  minutes  taken  for  the  required  number  of 

iterations. 

In  the  third  group  of  m  figures,  m  is  the  number  of  roots  of  the  original 
correlation  matrix  greater  than  unity.  The  m  figures  in  this  group  are  the  m 
largest  roots  of  the  matrix  for  the  specified  seeding  and  loss  function  parameters. 

12.2  The  Second  Line  or  Sequence  of  Lines 

The  second  line  or  sequence  of  lines  consists  of  four  numbers  to  a  line.  The 
numbers  all  have  to  do  with  the  simple  structure  transformation  described  in 
Chapter  9*  Since  data  sets  11  and  12  have  only  a  single  factor,  no  transformations 
were  required,  hence  no  lines  of  four  numbers  appear  for  these  data  sets.  For 
data  sets  1  through  9,  no  more  than  one  or  two  lines  are  given.  For  data  set  10, 

20  lines  of  four  numbers  each  are  given.  The  four  numbers  of  each  line  have  tie 

9 

following  interpretations: 

;  (l)  The  first  number  is  the  quantity  tr  (DA)  where  D  and  A  are  defined  in 

Eqs.  9*21  and  9*3^*  This  value  is  calculated  at  each  iteration  for  the  simple 
structure  transformation  matrix.  When  these  values  for  two  successive  iterations 
are  within  the  specified  tolerance  limit,  the  iterations  cease.  An  iteration  limit 
is  also  specified  in  the  computer  program  beyond  which  iterations  cease  even  though 
the  tolerance  limit  is  not  yet  reached. 

(2)  The  second  number  is  the  simple  structure  criterion  ¥  given  in  Eq.  9*70. 
The  maximum  value  this  criterion  can  attair  is  unity. 

(3)  The  third  number  is  the  number  of  sets  of  iterations  taken  to  calculate 
the  simple  structure  factor  loading  matrix  for  a  given  positive  integer  W  used  to 
calculate  F  in  Eq.  9*58.  This  integer  is  1  for  the  first  set  of  iterations.  If 
the  number  of  negative  factor  loadings  in  any  column  is  less  than  the  number  of 
factors,  the  computations  cease.  If  not,  the  integer  W  is  increased  by  1  and  a 

P 

second  set  of  iterations  for  the  simple  structure  factor  loading  matrix  occurs. 
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This  procedure  continues  until  at  least  one  column  of  the  simple  structure  factor 
loading  matrix  has  fewer  negatives  than  the  number  of  factors  or  columns  in  the 
matrix.  When  this  occurs,  the  simple  structure  matrix  fromlthe  preceding  set  of 
iterations  is  taken  as  the  final  simple  structure  factor  matrix,  except  for  W  *  1, 
in  which  case  the  corresponding  simple  structure  factor  loading  matrix  is  accepted. 

A  limit  is  put  on  the  number  of  successive  sets  of  iterations.  If  this  limit  is 
reached  before  the  number  of  negative  values  in  any  column  of  the  simple  structure 
factor  loading  matrix  is  less  than  the  number  of  factors,  the  successive  sets  of 
iterations  cease. 

(4)  This  number  is  the  integer  W.  It  indicates  the  number  of  sets  of  itera¬ 
tions  calculated  and  is  therefore  the  number  of  the  line  in  the  sequence  of  lines. 
The  integer  W  serves  as  the  argument  for  the  power  function  in  the  numerator  terms 
of  the  criterion  function  Y  given  by  Eq.  9-70  and  efined  in  more  detail  in  preced¬ 
ing  parts  of  Chapter  9. 

12.3  The  Factor  Loading  Lines 

A  set  of  n  lines,  where  n  is  the  number  of  variables,  is  given  for  each  of 
the  data  sets.  The  columns  in  this  set  of  lines  are  as  follows: 

(1)  The  first  column  gives  simply  the  line  numbers  which  indicate,  of  course, 
also  the  arbitrary  serial  numbers  of  the  variables  or  tests  described  in  Chapter  12. 

(2)  The  second  column  has  a  1  if  the  variable  retains  its  original  sign  and 
-1  if  its  sign  is  reversed  as  discussed  in  Chapter  6,  Section  4. 

(3)  The  third  column  gives  the  communalities  of  the  variables  as  calculated 
from  the  facto*-  loading  matrix  calculated  by  the  methods  of  Chapter  8. 

(4)  The  fourth  column  gives  the  specificities  corresponding  to  the  commun¬ 
alities  in  the  second  column.  The  sum  of  corresponding  elements  of  the  two  columns 
is  therefore  unity. 
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(5)  The  next  block  of  m  columns  gives  the  factor  loading  matrix  for  m  factors 
calculated  by  the  methods  of  Chapter  8. 

(6)  The  last  block  of  a  columns  for  data  sets  1  through  10  gives  the  simple 
structure  factor  loading  matrix  calculated  by  the  methods  of  Chapter  9*  For  data 
sets  11  and  12,  this  block  of  columns  is  omitted  since  only  one  factor  loading' 
vector  was  calculated  for  each  set. 

12.4  The  Last  Line 

The  last  line  for  data  sets  1  through  10  consists  of  m  figures,  where  m  is  the 
number  of  factors  solved  for.  Each  value  is  the  corresponding  element  of  the  A 
diagonal  matrix  defined  in  Eq.  9*34.  These  elements  are  the  ratios  whose  average 
is  given  by  Y  in  Eq.  9. TO.  It  is  the  average  of  these  ratios  which  is  the  second 
number  in  the  second  row  or  sequence  of  rows  described  in  Section  12.2  This  is 
the  simple  structure  criterion  we  seek  to  maximize.  The  maximum  value  any  one  of 
these  numbers  can  take  is  unity. 

This  final  row  of  figures  is  not  given  for  data  sets  10  and  11,  since  only 
one  factor  vector  was  calculated  f *r  each. 
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CHAPTER  13 


DISCUSSION  AND  CONCLUSIONS 

13*1  The  Six  Special  Cases 

Before  discussing  the  results  of  the  analyses  of  the  twelve  data  sets  it  may 
be  useful  to  relate  the  special  cases  of  the  loss  and  scaling  function  parameters 
to  traditional  factor  analysis  techniques.  For  convenience,  we  shall  construct 
six  five -letter  acronyms.  The  first  two  letters  will  indicate  the  scaling  parameter 
and  the  last  three  the  loss  parameter. 

If  the  scaling  parameter  is  p  =  0,  the  scaling  matrix  involves  only  the 
HEsiajal  variance  and  the  first  two  letters  of  the  acronym  for  this  case  will  be 
RE.  If  the  scaling  parameter  is  p  =  . 5,  the  scaling  matrix  involves  the  sum  of  the 
residual  and  estimated  or  common  variances  which  is  the  TOtal  variance,  and  the 
first  two  letters  of  the  acronym  for  this  case  will  be  TO.  If  the  scaling  parameter 
is  p  =  1,  the  scaling  matrix  involves  only  the  Estimated  or  common  variances  and  the 
first  two  letters  of  the  acronym  for  this  case  will  be  ES. 

If  the  loss  parameter  is  =  1,  the  loss  function  involves  only  the  residual 
COVariances  and  the  last  three  letters  of  the  acronym  for  this  case  will  be  COV. 

If  the  loss  parameter  is  »  0,  then  the  loss  function  involves  both  residual 
VAriances  and  Covariances,  equally  weighted,  and  the  last  three  letters  of  the 
acronym  for  this  case  will  be  VAC.  Thus  ve  shall  have  the  following  cases: 


p  *  0,  pw  *  1,  RECOV 

p  =  . 5,  PM  •  1.  TOCOV 

p  =  1.,  Py  =  1.  ESCOV 

p  -  0,  Pw  3  0  kevac 

p  *  .5,  py  3  0  TOVAC 

P  ■  1.,  Py  «=  0  ESVAC 


RECOV  factor  analysis  it  closely  associated  with  maximum  likelihood  factor 
analysis  developed  by  Levley  (19*40)  and  with  canonical  factor  analysis  developed 
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by  Rao  (1955)*  As  a  matter  of  fact,  the  equations  to  be  satisfied  by  RECOV  and 
the  latter  two  are  equivalent.  TOCOV  factor  analysis  is  closely  associated  with 
minres  factor  analysis  developed  by  Harmon  (1967).  ESCOV  factor  analysis  is 
similar  to  alpha  factor  analysis  developed  by  Kaiser  and  Caffrey  (I9S5).  R2YAC 
factor  analysis  has  been  discussed  by  Anderson  and  Rubin  (1956)  who  have  pointed  out 
fundamental  difficulties  with  the  model.  These  we  have  met  by  the  imposition  of 
somewhat  unsophisticated  ccnputatlmal  constraints.  TOVAC  factor  analysis  is  the 
same  as  what  many  investigators  call  principal  components  analysis.  Actually,  any 
of  the  methods  of  factor  analysis  which  is  a  special  case  of  the  loss  and  scaling 
parameters  we  have  discussed  may  be  regarded  as  a  principal  components  analysis 
of  a  real  symmetric  matrix,  whether  or  not  all  its  roots  are  non-negative.  ESVAC 
factor  analysis,  to  our  knowledge,  has  been  previously  discussed  only  by  the  author 
(Horst,  1965). 

13.2  Summary  of  Results 

We  recall  that  m  is  the  number  of  roots  of  the  correlation  matrix  greater 
than  unity  and  that  for  each  data  set  this  was  the  number  of  factors  solved  for. 

Vte  recall  also  that  the  criterion  4>  is  the  ratio  of  the  sum  of  squares  of  the  m 
largest  roots  of  a  matrix  with  specified  less  and  scaling  parameters  to  t:.e  tom 
of  squares  of  all  of  its  roota.  The  iterations  for  4  continued  until  the  absolute 
value  of  the  difference  between  tvo  successive  $'s  was  less  than  •  C'K'Cl  v:t>.  tr, 
upper  limit  of  100  iterations.  The  f  criterion  of  simple  structure  is  defined  in 
Chapter  9*  The  tolerance  limit  for  the  simple  structure  criterion  was  approximate) v 
the  same  as  for  The  actual  function  used  for  the  tolerance  .limit  is  discussed 
in  Chapter  12,  Section  2.  The  iteration  limit  was  60.  The  limit  on  the  number  cf 
sets  of  Iterations  was  20. 

Table  1  summarizes  results  for  the  twelve  data  sets. 

The  first  column  gives  simply  the  arbitrary  serial  numbers  of  the  data  se*\~. 
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The  second  column  gives  abbreviated  identification  of  the  data  sets.  The  third 
column  gives  the  number  of  factors  solved  for.  The  next  column  headed  "Grit." 
(Criterion)  gives  for  each  data  set  first  the  approximation  criterion  4  and  below 
it  the  simple  structure  criterion  f.  To  the  right  of  each  $  are  respectively  the 
number  of  iterations  and  the  actual  $  value  for  the  six  combinations  of  loss  and 
scaling  parameters.  The  first  six  columns  following  the  column  of  criterion 
symbols  give  the  data  for  ?w  =  1  and  the  next  six  columns  for  P^  *  0.  Within 
these  sets  of  six  columns,  the  first  pair  of  columns  gives  the  data  for  p  =  0  and 
the  next  two  for  p  -  .5  and  p  =  1  respectively. 

13.3  Rankings  by  Simple  Structure  Criterion 

Perhaps  one  of  the  most  important  queotions  to  be  answered  is  which  of  the 
six  methods  of  analysis  is  the  best  as  Judged  from  the  analysis  of  the  twelve  sets 
of  data.  One  overall  standard  might  be  based  on  the  simple  structure  criterion  T. 
Table  2  gives  for  each  of  the  first  nine  data  sets  the  rank  order  of  the  f  value 
for  each  of  the  six  combinations  of  loss  and  scaling  parameters.  The  first  three 
columns  of  rankings  are  for  P^  -  1  ar.d  the  last  3  for  Py  0  0.  Within  each  set 
of  three  columns,  the  first  gives  the  ranking  for  p  »  0,  and  the  next  two  for 
p  •  .5  and  p  *  1.  Rankings  for  wily  the  first  nine  of  the  twelve  data  sets  are 
given,  since  the  T’s  for  data  set  10  are  all  equal  and  no  T's  are  given  for  data 
sets  11  and  12. 

The  t  row  following  tne  row  for  data  set  )  gives  tbe  sum  of  the  rankings  for 
each  pair  of  less  and  scaling  parameters.  The  next  row  gives  the  rank  order  of 
the  sum  of  the  rankings  of  these  sums.  Ifce  next  row  H  gives  the  number  of  data 
sets  for  which  T  had  the  highest  ranking,  and  the  last  row  L  gives  the  number  of 
data  sets  for  which  T  had  the  lowest  rank. 

It  is  clear  from  the  last  three  rows  of  this  table  that  the  best  method  accord¬ 
ing  to  the  simple  structure  criterion  is  ESCOV  for  p  -  0  and  Py  •  1.  This  is  the 
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alpha  factor  analysis  model  of  Kaiser  and  Caffrey  (1965).  The  poorest  method 

is  KEVAC  for  p  =  0  and  P  =  0.  This  is  the  model  discussed  by  Anderson  and  Rubin 

w 

(I956).  In  view  of  the  problems  encountered  in  this  model,  it  is  not  surprising 
that  it  is  poorest  according  to  the  simple  structure  criterion. 

For  these  nine  data  sets,  RECOV  with  p  =  0,  P„  -  1  is  second  poorest  according 

w 

to  rankings  of  the  simple  structure  criterion.  This  is  the  method  closely  related 
to  the  maximum  likelihood  and  canonical  models  of  Lavley  (I9h0)  and  Rao  (1955) 
respectively. 

The  methods  in  second  and  third  place  respectively  are  TOCOV  with  p  =  .5 
and  ?w  =  1,  and  TOVAC  with  p  =  .5  and  =  0.  These  correspond  respectively  to  the 
minres  model  of  Harmon  (I967)  and  to  the  classical  principal  components  model. 

Obviously,  of  course,  the  procedure  we  have  used  for  the  comparative  evalua¬ 
tion  of  the  six  models  is  crude  and  is  based  on  a  very  limited  number  of  data  sets. 
Furthermore,  the  criterion  for  the  number  of  factors  is  arbitrary  and  other  criteria 
may  yield  different  results.  In  any  case,  it  is  quite  possible  that  for  any 
particular  data  set  one  may  wish  to  determine  the  loss  and  scaling  parameters 
p  and  Pw  so  as  to  maximize  the  criterion  Y.  Such  a  procedure  need  not  limit  the 
value  of  these  parameters  to  those  used  in  this  analysis. 

13.4  Ranks  by  Number  of  Iterations  Required 

A  further  procedure  for  a  relative  evaluation  of  the  six  models  may  be  based 
on  the  number  of  iterations  required  for  0  to  stabilize  to  the  specified  tolerance 
limit.  Table  3  provides  an  analysis  similar  to  that  of  Table  2.  Here,  however, 
the  rankings  are  for  the  number  of  iterations  required  for  0  given  in  Table  1  and 
all  12  data  sets  sure  included.  In  this  ranking  we  exclude  the  2-2  column  since 
this  is  the  principal  axis  method  and  except  for  peculiarities  of  the  computer 
program,  no  iterations  would  be  required.  In  any  case,  according  to  the  last  three 
ro ws  of  Table  3>  ESVAC  with  p  ■  1,  *  0  is  best  and  ESCOV  with  p  =  1,  P^  =  1  is 
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second  best.  This  latter  is  the  model  which  came  out  best  in  Table  2  based  on 
the  simple  structure  criterion  Y.  The  poorest  is  HEC0\  with  p  =  0  and  =  1, 
which  came  out  second  poorest  in  Table  2. 

Again  it  is  obvious  that  the  rating  procedure  is  crude  and  based  on  what  some 
may  regard  as  a  questionable  criterion  for  the  number  of  factors.  However,  a 
cursory  examination  of  Table  3  shows  that  in  general  the  number  of  iterations 
required  for  0  to  stabilize  tends  to  be  substantially  greater  for  RECOV  than  for 
the  other  models.  It  is  believed  that  this  tendency  would  persist  even  with  other 
defensible  criteria  for  the  number  of  factors. 

In  view  of  the  marked  increase  in  iterations  required  for  this  model  over 
those  required  for  other  models  and  at  least  some  persuasive  indications  of  poor 
simple  structure  potentiality,  one  might  question  whether  this  model  is  to  be 
generally  recommended.  Since  it  is  closely  related  to  the  maximum  likelihood  and 
canonical  models,  one  might  also  question  whether  the  great  interest  and  effort 
accorded  these  models  in  the  past  is  completely  justified. 

13*5  The  Simple  Structure  Factors 

A  more  detailed  examination  cf  one  simple  stiucvure  factor  matrices  given  at 
the  right  of  the  tables  in  Chapter  12  for  data  sets  1  through  10  might  be  of 
interest.  No  simple  structure  matrices  were  calculated  for  data  sets  11  and  12 
since  only  one  factor  was  obtained.  The  largest  element  in  each  row  is  underline-*- 
for  all  the  simple  structure  matrices  in  data  sets  1  through  10.  For  each  data 
set  there  are  six  of  these,  one  for  each  combination  of  the  loss  and  scaling 
parameters.  (See  Chapter  12  for  detailed  description.)  For  the  first  nine  data 
sets  these  simple  structure  matrices  may  be  compared  with  those  obtained  by  other 
investigators  referred  to  in  Chapter  11.  Brief  comments  on  the  data  sets  might  be 
of  interest  as  follows: 
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1.  Primary  Mental  Abilities.  For  all  six  models,  the  simple  structure 
matrices  are  sharp  and  agree  well  with  the  results  of  Thurstone  and  Thurstone 
(1941). 

2.  Twenty -four  psychological  tests.  According  to  the  Y  criterion,  ESCOV 
with  p  =  1,  P =  1  gives  the  best  simple  structure.  This  model  corresponds  to 

Kaiser's  (1965)  alpha  model.  A  number  of  simple  structure  solutions  are  given  by 
Harmon  (1967)  for  this  data  set.  However,  his  various  solutions  involve  only  four 
factors  whereas  ours  has  five.  But  all  of  his  simple  structure  solutions  give 
results  easily  recognized  as  similar.  Our  simple  structure  for  ESCOV  gives  results 
similar  to  his  for  three  of  the  factors.  However,  his  solutions  assign  variables 
14  through  19  essentially  to  a  single  factor,  whereas  ESCOV  splits  them  into  two 
factors,  the  first  three  going  to  one  and  the  last  three  to  another  as  follows : 

14.  Word  recognition  17*  Object  -  number 

15.  Number  recognition  18.  Number  -  figure 

16.  Figure  recognition  19.  Figure  -  word 

Without  referring  in  detail  to  a  description  of  the  original  tests,  it  is  not 
surprising  .hat  recognition  of  various  types  of  visual  stimuli  should  have  a  factor 
in  common,  and  ability  to  associate  pairs  involving  two  different  types  of  stimuli 
should  have  another  factor  in  common.  It  is  quite  possible,  of  course,  that  if 
Harmon  had  included  a  fifth  factor  in  his  analysis  he  also  would  have  found  the 
same  factor  differentiation  between  these  two  triplets  of  tests.  It  is  interesting 
to  note  that  our  TOCOV  solution  which  corresponds  to  Harmon's  rainres  does  not 
appear  to  give  as  sharp  a  simple  structure  for  five  factors  as  his  does  for  the 
four  factors  on  which  he  uses  a  direct  oblimin  solution. 

3.  Thirty-three  variable  speed  study.  Of  the  six  models,  ESCOV  again  gives 
the  highest  simple  structure  criterion  for  the  data  from  Lord  (1956).  Lord  solved 
for  ten  factors  by  a  modification  of  the  maximum  likelihood  method.  Our  criterion 
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for  the  number  of  factors  yielded  only  six  factors.  Lord's  rotations  were  carried 
out  by  subjective  non -analytical  procedures  so  that  his  simple  structure  matrix 
is  not  comparable  to  our  ESCOV  with  P„  =  1  and  p  =  1.  However,  referring  to 
Section  11.3,  the  verbal,  spatial,  and  number  factors  come  out  clearly  as  they  do 
in  Lord's  analysis.  In  addition,  a  factor  common  to  the  two  number  speed  tests  (23) 
and  (24)  and  Lord's  reference  tests  for  perceptual  speed  (25),  (26),  and  (27) 
appears  in  our  analysis.  It  is  also  interesting  that  in  our  analysis  college  grades 
tend  to  split,  with  English  (28)  going  to  the  verbal  factor,  Engineering  Drawing 
(30)  to  the  spatial  factor,  and  Foreign  Language  (29),  Chemistry  (31),  and  Mathe¬ 
matics  (32)  predominately  to  a  factor  thrt  we  might  characterize  as  a  facility  with 
symbolic  systems.  A  factor  which  Lcrd  failed  to  find  has  only  a  single  high  load¬ 
ing  on  Conduct  (33)  and  small  positive  loadings  for  the  five  grade  variables  (29) 
through  (32)  This  may  be  a  conformity  factor  considering  the  data  are  based  on 
students  in  a  military  academy. 

4.  Thurstone's  twenty -variable  box  problem.  From  this  classical  set  of  data 

there  is  little  to  choose  among  the  simple  structure  factor  matrices  for  the  six 
models.  As  would  b®  ''xpec+od,  three  factors  were  obtained.  With  the  exception  of 
RECOV  for  p  =0,  =  1,  the  ¥  values  do  not  differ  by  more  than  .001.  For  HECOV 

the  value  is  only  .003  less  than  the  highest  value  of  .646  for  TOCOV,  ESCOV,  TOVAC, 
and  ESVAC.  The  simple  structure  is  clear  for  all  cases  with  the  X,  X,  and  Z 
dimension  variables  coming  out  with  only  a  single  _arge  loading  for  a  factor  and 
the  functions  of  these  variables  indicated  in  Section  11.4  having  the  loadings 
that  would  be  expected. 

5.  Eight -variable  body  type  measures.  All  six  models  for  this  data  set, 
yielding  only  two  factors,  give  very  clear  simple  structure  as  has  been  found  by 
other  investigators. 


— *<***> **»**■  , i,  nrwwwM  —  -“**« 
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6.  Twelve -variable  anthropometric  measures.  In  an  early  analysis  of  these 
data,  Thurstone  (1946)  found  four  simple  structure  factors  by  subjective  graphical 
methods  whereas  our  criterion  gave  only  three  factors.  For  all  models  except 
RECOV  and  REVAC  whose  loss  parameter  p  is  0,  the  Y  values  are  in  the  high  .80's 
with  ESCOV  again  being  highest  with  Y  =  .880.  The  simple  structures  for  all  six 
models  are  clear  cut.  Thurst one's  B  and  D  factors  have  both  loadings  of  .45  or 
more  on  the  three  variables  stature  (l),  span  (5),  and  hand  length  (11).  In 
addition,  his  factor  B  has  high  loadings  on  sitting  height  (2),  and  his  factor  D 
has  a  high  loading  on  hand  breadth  (12).  Our  first  factor  tends  to  collapse  these 
two  factors,  while  our  second  and  thifd "factors  are  easily  recognized  as_Thurstone 
factors  C  and  A  respectively. 

7*  Fifteen  variables  from  Hemmerle.  Hemmerle  (1965)  does  not  give  the  source 
of  this  matrix  nor  does  he  identify  the  variables.  Since  he  was  concerned  primarily 
with  a  computational  procedure  for  the  maximum  likelihood  factor  model,  he  did  not 
attempt  a  simple  structure  transformation.  He  extracted  eight  factors  but  does  net 
state  his  criterion.  Cur  criterion  got  only  five  factors.  These  data  were  included 
in  our  study  primarily  because  Joreskog  (1967)  as  well  as  Horst  (1968b)  had  also 
worked  extensively  with  them.  Both  of  us  had  found  the  data  to  behave  peculiarly 
and  our  earlier  results  for  eight  variables  were  markedly  different  for  the  maxi¬ 
mum  likelihood  methods  of  Joreskog  and  our  corresponding  RECOV  model.  Even  in  the 
present  study  it  is  the  only  data  set  that  has  its  highest  Y  value  (.596)  for  REVAC 
which,  from  a  theoretical  point  of  view,  is  the  poorest  of  the  six  models.  Since 
nothing  is  known  about  the  identity  of  the  variables,  nothing  of  substantive 
interest  can  be  said  about  their  simple  structure  factor  loadings. 

8.  Seventeen -variable  data  from  Bechtold- -Sample  1.  For  this  data  set  our 
criterion  yielded  only  five  factors  whereas  Bechtold  (1961)  had  deliberately 
attempted  to  represent  six  factors  in  his  battery,  as  indicated  in  Section  11.8. 
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In  general,  his  V,  W,  S,  and  N  factors  cone  cut  clearly  in  all  bIa  models.  The 
memory  factor  for  variables  1  and  2  fails  to  come  out  clearly  in  EECOV  which  has 
the  highest  Y  value,  although  the  reasoning  factor  R  for  variables  15,  1 6,  and  IT 
comes  out  clearly  in  this  model.  The  only  other  model  for  which  R  comes  out 
clearly  is  the  suspect  REVAC  which  also  fails  on  the  M.  It  is  q.uite  probable  that 
our  criterion  for  number  of  factors  was  too  low  for  this  data  set. 

9*  Seventeen -variable  data  from  Bechtold- -Sample  2.  Since  the  tests  in  this 
data  set  were  the  same  as  for  data  set  8  and  the  sample  was  presumably  comparable, 
the  results  should  be  substantially  the  same  for  the  s'  -)le  structure  matrices. 

As  in  the  previous  set,  the  criterion  yielded  five  factors.  For  all  six  models, 
the  V,  W,  S,  and  N  factors  are  clearly  defined.  As  in  the  previous  set,  the  M 
factor  for  variables  1  and  2  appears  in  all  models  except  RECOV  a'  1  REVAC  but  some¬ 
what  less  clearly  in  TOCOV.  The  R  factor  for  variables  15,  16,  and  17  appears  most 
clearly  as  a  distinct  factor  in  RECOV  and  TOCOV.  Here  again,  it  is  highly  probable 
that  our  criterion  for  number  of  factors  was  too  restrictive. 

10.  Nine -variable  synthetic  data.  The  origin  or  source  of  this  data  set  is 
described  in  detail  in  Chapter  11,  Section  10.  The  correlation  matrix  was  con¬ 
structed  so  that  the  simple  structure  factor  loading  matrix  would  have  three 
factors  and,  to  three -decimal  accuracy,  this  matrix  would  be  as  follows: 


I 

II 

III 

1 

000 

.383 

.921* 

2 

.000 

.707 

.707 

3 

.000 

.92I* 

.383 

k 

.383 

924 

.000 

5 

'  .TOT 

.707 

.000 

6 

.921* 

.383 

.000 

T 

,0°.h 

.000 

.383 

8 

.707 

.000 

.707 

9 

.383 

.000 

.921* 

It  can  be  shown  that  if  one  allocates’  three  points  on  each  of  the  arcs  of  a 
right  spherical  triangle  as  indicated  in  Chapter  11,  Section  10,  then  to  three 
decimal  places  the  cosines  of  the  angles  of  eacn  of  these  nine  points  with  each  of 
the  three  vertices  of  the  spherical  triangle  will  be  as  shown  in  the  matrix  above. 
For  some  decades  we  have  been  trying  to  find  an  analytical  procedure  for  recovering 
this  matrix  from  the  correlation  matrix  of  these  points.  To  our  knowledge,  none  of 
the  analytical  methods  previously  available  will  accomplish  this  recovery.  Refer - 
er.  a  to  Chapter  12  for  the  results  from  this  data  set  shows  that  the  simple  struc¬ 
ture  factor  matrices  for  all  six  models  differ  at  most  from  the  above  matrix  by 
.001.  For  all  six  models,  the  sets  of  simple  structure  iterations  went  to  the 
prespecified  limit  of  twenty.  For  each  model  the  number  of  iterations  for  the 
first  set  went  to  the  prespecified  limit  of  sixty.  Thereafter,  however,  the  number 
of  iterations  required  for  the  successive  sets  diminished  rapidly  to  two  or  three. 
It  is  quite  probable  that  if  no  limit  were  placed  on  the  number  of  sets  of  itera¬ 
tions,  the  original  simple  structure  matrix  could  be  recovered  to  any  desired  degr< 
of  accuracy. 

13.6  Improper  Solutions 

A  factor  analysis  model  whose  loss  function  involves  only  the  residual  covari¬ 
ances  may  yield  communalitles  for  some  tests  which  exceed  unity.  Such  a  result  ic 
known  as  a  Heywood  case.  These  factor  analysis  models  are  our  RECOV,  TOCOV,  and 
ESCOV  which  correspond  respectively  to  Lawley's  (19^0)  maximum  likelihood,  Harmon's 
(1967)  minres,  and  Kaiser’s  (1965)  alpha.  It  ^ 3  of  interest  to  note  that  for  none 
of  our  data  sets  does  our  method  of  computation  involving  real  data  give  commun- 
alities  as  high  as  the  constrained  values  of  .9995  for  any  of  these  three  cases. 
This  is  not  true  for  Harmon’s  (19&7)  minres  on  data  set  2,  where  he  obtains  his 
maximum  constrained  cownunality  of  unity  for  variable  19. 
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When  an  unconstrained  solution  yields  camnunalities  greater  than  unity,  this 
is  sometimes  called  an  improper  solution.  In  the  case  of  maximum  likelihood  solu¬ 
tions,  Joreskog  (196?)  observes  that,  "Experience  varifies  that  improper  solutions 
are  found  more  often  than  is  usually  expected. "  Our  RECOV  computational  algorithms 
have  been  applied  to  some  of  the  data  sets  on  which  Joreskog  has  applied  his  compu¬ 
tational  algorithms  for  the  corresponding  maximum  likelihood  method.  In  general, 
where  we  have  used  the  same  number  of  factors,  neither  of  us  has  encountered  an 
unconstrained  improper  solution.  However,  for  the  case  of  data  set  7  from  Hemmerle 
for  eight  factors,  Joreskog' s  (1967)  procedure  found  it  necessary  to  constrain 
variables  7  and  15  whereas  our  (Horst,  1968b)  procedure  up  to  10,000  iterations 
found  no  improper  communalities.  A  highly  accelerated  modification  of  our  pro¬ 
cedure  did  not  require  the  constraining  supplements  in  the  computational  procedures. 
Nevertheless,  it  is  not  improbable  that  for  a  completely  adequate  factor  analysis 
system,  the  occurrence  of  improper  communalities  would  signal  either  the  use  of  in¬ 
appropriate  loss  or  scaling  parameters,  an  inappropriate  criterion  for  the  number 
of  factors,  or  some  combination  of  these  J 
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COMPUTER  PROGRAMS 


The  Fortran  IV  computer  programs  for  carrying  out  the  analyses  of  the  previous 
chapters  consist  of  a  main  program  (MAIN)  and  overlay  subroutine  subprograms  called 
by  the  main  program.  The  overlay  subroutines  are  called  SYMI,  JACS3,  JACS,  RARE, 
SIMP;  and  IUPLI . 

14.1  MAIN 

The  main  program  provides  parameter  values  required  for  the  computations,  an 
outer  loop  for  the  data  sets,  an  intermediate  loop  for  the  loss  parameter,  an 
inner  loop  for  the  scaling  parameter,  and  a  call  to  the  output  overlay  subroutine 
BUPLI. 


The  parameters.  It  is  standard  practice  to  read  in  parameters  from  cards  along 
with  the  data  cards.  Particularly  is  this  true  if  the  program  deck  is  a  binary 
deck.  It  is  our  opinion  that  binary  program  decks  are  essentially  obsolete,  especi¬ 
ally  with  the  rapid  compilers  currently  available.  It  is  usually  desirable  to  have 
the  source  program  immediately  available  with  the  output  of  a  given  computer  run, 
together  with  all  the  program  parameters  and  option  codes  that  were  used  in  the 
computer  run.  We  have  been  repeatedly  frustrated  in  attempting  to  assist  laymen  in 
the  interpretation  of  their  computer  output  by  the  fact  that  they  used  binary  pro¬ 
gram  decks  and  therefore  could  provide  no  information  about  fhe  program  parameters, 
option  codes,  and  the  computing  algorithms  utilized. 

If  a  source  program  deck  such  as  Fortran  IV  is  used,  it  is  possible  to  read 
in  program  parameters  and  option  code  cards  as  data  and  these  cards  can  be  varied 
to  suit  the  requirements  of  the  investigator  and  his  data.  However,  it  may  be  con¬ 
venient  in  research  with  various  data  analysis  models  to  provide  some  of  this 
information  in  program  statements  so  that  they  may  be  readily  found  at  the  beginnin> 
of  the  program  listing.  Some  of  these  values  are  given  at  the  beginning  of  MAIN. 

A  number  of  them  are  repeated  vith  different  numberical  values.  The  last  time  the 
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replacement  statement  appears  for  a  parameter  variable  is  of  course  the  value  it 
takes  in  the  program.  It  has  been  found  convenient  for  research  purposes  to  pro¬ 
vide  a  number  of  values  which  may  be  changed  merely  by  changing  the  position  of 
the  statement. 

The  parameters  area  as  follows: 

P  tolerance  limit 

LIB  beginning  indexing  parameter  for  loss  parameter 

LIE  ending  indexing  parameter  for  loss  parameter 

LB  beginning  indexing  parameter  for  scaling  parameter 

IE  ending  indexing  parameter  for  scaling  parameter 

NL  iteration  limit  for  simple  structure  iterations 

NF1  ending  indexing  parameter  fear  row  scaling  option  of  factor  matrix 

KKL  iteration  limit  for  principal  axis  solution 

EE  Tolerance  limit  for  specificity  variance 

The  outer  loop.  This  is  the  loop  with  the  index  LLL  and  the  indexing  para¬ 
meters  1,  NP.  This  loop  controls  the  number  of  data  sets  processed  in  a  given 
run.  The  data  for  each  set  consists  of  the  number  of  variables  in  the  set,  the 
format  of  the  correlation  matrix,  and  the  correlation  matrix,  itself.  The  loop 
calls  the  subroutine  SYMI  which  provides  initial  estimates  of  the  residual  variance 
and  JACS3  which  determines  the  number' of  factors. 
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The  intermediate  loop.  This  is  the  loop  with  the  index  LLI  and  indexing 

parameters  LIE,  LIE.  It  calculates  the  loss  parameter  P .  In  this  program  the 

w 

calculation  of  only  P„  =  1  and  P„  =  0  are  provided  for  but  any  desired  intermediate 

W  w 

values  could  be  provided  with  slight  modification. 

The  inner  loop.  This  loop  has  the  index  LL  and  the  indexing  parameters  LB, 

LE.  It  ■'■rites  the  parameters  ILL,  N,  LI,  LLI,  LL,  and  NF1  on  scratch  tape.  It 
calculates  functions  of  the  scaling  parameter  p  which  are  used  in  the  calculation 
of  the  scaling  matrix.  In  this  program  only  three  loss  function  parameters  are 
provided  for.  These  are  p  =  0,  p  =  .5,  and  p  =  1.  However,  as  in  the  case  of  the 
loss  parameter,  any  desired  intermediate  values  could  be  provided  with  slight  modi¬ 
fication. 

This  loop  also  calls  JACS  which  calculates  a  first  approximation  to  the  basic 
structure  factor  matrix  and  RARE  which  calculates  iteratively  the  simple  descaled 
matrix  from  the  basic  structure  factor  loading  matrix.  If,  as  is  usually  the  case, 
the  number  of  factors  exceeds  1,  this  loop  also  calls  SIMP  which  calculates  the 
simple  structure  factor  loading  matrix. 

KJPLI.  This  subroutine  is  outside  the  outer  loop  of  MAIN.  It  reads  from 
scratch  tape  the  data  that  is  to  be  printed  and  writes  it  on  BCD  tape  in  the  format 
in  which  the  data  in  Chapter  12  are  given. 

14.2  SYMI 

This  overlay  subroutine  reads  in  the  data,  calculates  the  inverse  of  the 
correlation  matrix,  and  then  calculates  the  first  approximation  to  the  residual 
variances. 

Data  input.  A  single  card  giving  the  number  of  variables  N  is  read  with  forma'1 
(l4).  The  program  has  been  dimensioned  for  up  to  80  variables.  It  could  probably 
be  extended  to  85  and,  with  some  rewriting,  to  90.  An  A-format  card  giving  the 
format  of  the  correlation  matrix  is  read.  The  correlation  matrix  is  read.  As  the 
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program  is  written,  each  row  of  the  correlation  matrix  must  begin  on  a  new  card. 

The  program  assumes  that  at  least  the  infra -diagonal  elements  are  given.  The 
supra -diagonal  elements  must  either  be  given  or  treated  as  zero.  In  either  case 
the  program  then  writes  the  supra -diagonal  elements  ana  enters  unity  in  the 
diagonals.  Then  the  correlation  matrix  is  stored  on  scratch  tape. 

Matrix  inversion.  The  program  calls  a  regular  subroutine  SYMIN  to  invert 
the  corre’ it ion  matrix.  If  SYMIN  finds  that  the  correlation  matrix  is  not  basic 
or  positive  definite,  it  returns  control  to  SYMI,  the  overlay  subroutine,  which 
shrinks  the  off diagonal  elements  by  a  factor  of  .9.  This  factor  is  arbitrary.  It 
can  be  6hown  that  if  the  off diagonal  elements  on  any  correlation  matrix  are  multi¬ 
plied  by  a  positive  value  less  than  unity,  the  resulting  matrix  will  be  basic  and 
hence  have  a,  regular  inverse.  SYMIN  is  again  called  to  invert  the  modified  corre¬ 
lation  matrix. 

Residual  variance  approximation.  The  reciprocals  of  the  diagonal  elements 
of  the  inverse  of  the  correlation  matrix  are  calculated.  These  provide  the  first 
approximation  to  the  residual  variances.  They  are  written  on  scratch  tape. 

14.3  JACS3 

This  subroutine  reads  the  correlation  matrix  from  scratch  tape  on  which  it 
was  written  by  SYMI.  It  calculates,  by  an  adaptation  of  the  Jacobi  method,  all 
the  roots  of  the  correlation  matrix  in  order  of  magnitude  which  are  greater  than 
unity.  It  transmits  the  number  of  these  roots  to  common  core  storage  as  the  number 
of  factors. 

14.4  JACS 

This  overlay  subroutine  reads  the  necessary  data  from  scratch  tape.  It  then 
calculates  the  first  approximation  to  the  modified  correlation  matrix  with  specified 
loss  and  scaling  parameters,  and  the  first  approximation  to  the  basic  structure 
matrix. 
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The  data.  The  correlation  matrix  is  read  from  scratch  tape  on  which  it  was 
written.  Without  rewinding,  the  first  approximation  to  the  residual  variance 
vector  is  read  from  the  same  tape. 

The  scaling  matrix.  The  communality  variance  is  calculated.  The  scaling 
matrix  is  calculated  as  a  function  of  the  communality  and  residual  variance  vectors 
and  the  scaling  parameter. 

The  modified  correlation  matrix.  The  correlation  matrix  is  scaled  by  the 
scaling  matrix.  The  diagonals  of  the  resulting  matrix  are  adjusted  according  to 
the  current  loss  parameter. 

The  basic  structure  matrix.  Subroutine  JACSIM  is  called.  This  subroutine 
calculates  the  first  m  principal  component  or  basic  structure  vectors  of  the  modi¬ 
fied  correlation  matrix  where  m  is  the  r umber  of  roots  greater  than  unity  found 
in  Section  14.3.  This  is  the  first  approximation  to  the  principal  axis  matrix  for  .. 
specified  loss  and  scaling  parameters.  The  principal  axis  matrix  is  written  on 
scratch  tape. 

14.5  RARE 

This  overlay  subroutine  calculates  the  descaled  principal  axis  matrix  for 
current  loss  and  scaling  parameters.  It  reads  the  necessary  data  from  scratch 
tape.  It  then  calculates  a  first  approximation  to  a  descaled  principal  axis  mati.'.:: 
and  a  second  approximation  to  the  scaling  matrix.  It  calculates  iteratively  the 
descaled  principal  axis  matrix  for  predetermined  loss  and  scaling  parameters.  It 
writes  output  data  on  scratch  tape.  Next  it  effects  row  sign  reversals  if  needuu. 
Finally,  it  transfers  daua  to  scratch  tape. 

Input  data.  The  correlation  matrix  is  read  from  scratch  tape.  The  first 
approximation  to  the  principal  axis  factor  matrix  is  read  frcm  another  scratch  tape. 

First  descaled  principal  axis  matrix.  The  descaling  matrix  for  the  first  ap¬ 
proximation  to  the  factor  loading  matrix  is  calculated.  The  first  approximation 
to  the  descaled  factor  loading  matrix  is  calculated. 
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The  scaling  matrix.  Second  approximations  to  the  commonality  and  residual 
variance  vectors  are  calculated.  A  second  approximation  to  the  scaling  vector  is 
calculated  from  these  two  vectors  and  functions  of  the  scaling  parameter. 

Successive  descaled  principal  axis  matrices.  A  loop  with  index  KKK  and  index¬ 
ing  parameters  1,  KKL  is  set  up  to  call  iteratively  subroutine  RARED.  This  sub¬ 
routine  calculates  successive  approximations  to  the  descaled  principal  axis  matrix  f 
for  the  loss  and  scaling  parameters  determined  within  the  inner  loop  of  MAIN.  The 
computations  are  carried  out  by  the  algorithms  indicated  in  Chapter  8.  The  sub¬ 
routine  includes  a  constraint  to  keep  the  residual  variances  positive.  It  also 
calculates  the  criterion  0  of  Chapter  12  and  the  difference  between  two  successive 
0’s  as  a  convergence  tolerance. 

Output  data.  The  final  $  value,  the  number  of  iterations  taken,  and  the  total 
time  in  seconds  are  written  on  the  scratch  tape  which  will  s'hsequently  be  read 
back  for  output.  The  first  m  roots  of  the  final  modified  correlation  matrix  are 
also  written  on  this  tape. 

Sign  reversals.  The  first  element  in  each  row  of  the  final  descaled  principal 
axis  matrix  is  checked  for  sign.  Sign  reversals  by  row  are  made  where  necessary. 

Transfers  of  output  data  to  tape.  The  final  descaled  principal  axis  matrix, 
together  with  the  sign  vector  and  the  final  communality  and  residual  variance 
vectors,  are  written  on  the  scratch  tape  for  output  data.  The  final  descaled 
principal  axis  matrix  is  also  written  on  another  scratch  tape  to  be  read  sube 
sequently  for  further  operations. 

14.6  SIMP 

This  overlay  subroutine  calculates  the  simple  structure  factor  loading  matrix 
It  provides  parameter  values  and  options  of  row  scaling  for  the  descaled  principal 
axis  matrix.  It  has  a  major  outer  and  an  inner  iteration  loop  for  calculating  the 
simple  structure  matrix.  It  writes  the  output  data  on  the  output  scratch  tape. 
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Parameter  values .  The  parameter  LLE  gives  the  limit  on  the  number  of  sets 
of  iterations.  The  parameter  ML  is  used  in  calculating  the  F  value  in  Eq.  9. 7* 

Rov  scaling  option.  The  program  normalizes  the  principal,  axis  matrix  by  rows 
before  beginning  the  simple  structure  iterations.  The  final  simple  structure 
matrix  is  denormalized  before  being  written  on  output  tape.  This  solution  is 
given  by  using  the  parameter  NF1  =  1  in  MAIN.  This  parameter  serves  as  the  end 
indexing  parameter  in  SIMP  for  the  DO  index  LLL.  If  in  addition  to  this  solution 
it  is  also  desired  to  have  a  solution  without  first  normalizing  by  rows,  the 
parameter  NF1  =  2  is  used  in  MAIN.  The  program  does  not  provide  for  just  the  non- 
normalized  solution  but  with  slight  modification  it  can  be  made  to  do  so. 

Outer  iteration  loop.  The  major  outer  iteration  loop  has  the  index  LLL4  with 
indexing  parameters  LLB,  LLE.  This  loop  p  yvides  for  successive  sets  of  iterations 
where  the  exponent  F  decreases  with  each  succeeding  set.  The  value  F  is  calculated 
as  a  function  of  the  index  LLL4  and  the  parameter  ML.  For  each  iteration  set,  this 
loop  also  writes  on  the  output  scratch  tape  the  tolerance  criterion,  the  simple 
structure  criterion,  the  number  of  iterations,  and  the  number  of  the  iteration  set. 
For  each  set  of  iterations,  this  loop  determines  whether  any  vector  of  simple 
structure  factor  loadings  has  less  than  m  negative  values.  If  so,  no  further  set 
of  iterations  is  calculated. 

Inner  iteration  loop.  This  loop  has  the  index  LL  with  the  indexing  paramete.rr 
1,  NL.  It  calculates  iteratively  the  transform'  ;ion  matrix  and  the  simple  structure 
factor  matrix  by  means  of  the  algorithms  given  in  Eqs.  9.57  through  9. 71.  For  each 
iteration  it  cells  the  subroutine  SYMI3  which  calculates  the  inverse  of  a  positive 
definite  symmetric  matrix.  Within  this  loop  also  is  calculated  the  criterion 
value.  If  two  successive  values  are  within  the  tolerance  limit,  the  iterations 


are  terminated. 
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The  final  simple  structure  matrix.  After  the  successive  sets  of  iterations 
are  terminated,  the  program  recognizes  as  the  simple  structure  factor  matrix  the 
one  calculated  in  the  next  to  last  set  of  iterations,  unless  only  one  set  was 
calculated.  In  the  latter  case,  the  matrix  calculated  in  the  single  set  of  itera¬ 
tions  is  recognized  as  the  simple  structure  matrix.  For  the  calculations  beginning 
with  a  row  normalization  of  the  principal  axis  factor  matrix,  the  final  simple 
structure  matrix  is  denormalized  by  rows.  In  either  case,  the  final  simple  struc¬ 
ture  matrix  is  written  on  the  output  scratch  tape.  The  vector  of  Y  criterion  values 
for  each  simple  stricture  factor  vector  is  also  written  on  the  output  scratch  tape. 
14.7  DUPLI 

This  overlay  subroutine  reads  the  data  on  the  output  scratch  tape  and  writes 
it  on  BCD  output  tape  according  to  the  format  of  the  data  in  Chapter  12.  The  sub¬ 
routine  has  an  outer  loop  with  index  LS  and  indexing  parameters  1,  NS  so  that  NS 
copies  of  the  output  will  be  printed. 
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MA  IN 

DIMENSION  J  f I  SO  J 

COMMON  P,NL,N,NF,LtKl  , KKL ,KK2L , NA ,F 1 » F F , KK3L ,HHf KKK 
*,NC»FF1,FF?,LL  1 
*tPD,OOfPW 
♦  .NR 
* »  T  I M 
*»LT 

* » L  L  A »  N°  t  L IP»L IEfLP  »LF»J! 

P  =  .00001 
NP=7 
NP  =  B 
NP=S 
NP=3 
NP=  1 
NP  =  4 
NP=? 

N°=  1  ? 

L  R=  3 
l  !«  =  ? 

LB  =  ? 

I  F=  ? 

L  0=  1 
L  IF  =  1 
l  I  H  =  1 
L  lc=^ 
l.  F  =  1 
I  F  =  3 
NL  =  ">0 
NL=  100 
NL  =60 
NR  =’ 

MFJ  =1 
KKL  =10 
KKL  -^O 
KKL  *  1 00 
Ff> .0 04 
Ft*  ,000r' 

MFWlMO  P, 

L  7*0 

WP IT-1 A|°rrU 

cr»Q  FORMAT!  IH1  1 

on  ppq  LLL  =1 ,NP 

call  $VM J 

CALI  JACS1 

00  «*?  U i-lth.l  IF 
PW*  “*-1  L  » 
on  ppi  t  l  *1  Of  L  r 

Wp  I  TF  I  l.'t  °0  7)LLLtNtLffi  I  ItlLfNF] 

90?  EOriMATH  I*  » 

WOITMCt  OQ>» 
cap  r OR MAl I / / 1 
l  7*1 7 ♦  1 

VfPITrjHM  I  LtN.LIfLlI  fLlfNFl 

00. T-t i 
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pr>=°n/2. 

Q0=  i.-po 
PQ=1.-'>.*P0*QD 
po=ph/pq 
00=QO/PQ 
C ML  JACS 
CALL  RAPT 
WR  ITT  ?) 

9P3  FORMAT!///) 

I  F( L I- 1 )RR0,R31 » 8  80 
880  CONTINUF 
CALL  SIMP 
J II L7  )=LL A 
PR l  CONTTNOF 
88?  CDNT  f  MUF 

WRITE  I  *,«<)<;) 

RRR  CONTI, NUF 
REWIND  « 

WR ITT! 6tT)Q| 

CALL  OUPLI 

STOP 

CN0 
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SUBROUTINE  SYMI 

DIMENSION  R I  80,801, Y( 80,80 1 jA I  80) 

*, DEI  80 1 
*»  EM  I  121 

COMMON  P,NL,N, NF,  L , KL , KKL ,KK2L ,NA , E 1 , E E , KK3L , HH, KKK 
*, NC»  FF1 ,FF2,LL  l 
*,PO,yD,PW 
* ,  NF  ! 

*,TIM 

*,l  I 

*,  L  L  A,  NP  ,L  IO,LIF,LR,LE,JI 
REWIND  2 
REWIND  3 
REWIND  4 
REAC! 5.992IN 
99?  FORMAT!  14  ) 

READ!  8,991  MFMIII  ,1  =1,  121 
991  FORMAT! 12A6I 
DO  50?  I  = 1 »N 

RFA015»FM||PI I,JI,J=1,NI 
502  CONTINUE 
DO  4  I  5 1 , N 
DO  2  J  *  I ,  N 
R!  I  ,J  »=R! J,  1 1 
2  CONTINUE 
R!  I, I  IM. 

4  CONTINUE 
no  1 7  I  *  1  ,  N 
W  R I T  E ( 2  I! R  <  I,  ,NI 

l?  continue 
I  s*o 

CALL  SYMIN!R,N, IS  I 
IF!  IS»40,44,6U 
60  CONTINUE 
REWIND  2 
00  6  1  I  •  l ,  M 
or-A.H  MIR!  I ,  J  I ,  J*  l ,  N I 

41  continue 

0(1  4'  1*1, N 
DO  6?  J  *  I  »  N 
«II,J1*P| I  ,  J  I  •  •  9 
4?  CONTINUE 
AlltlH. 

AT  CONTINUF 

call  syminir,n, m 

44  CONTIflUF 

DO  501  !■  1  ,N 
Of!  II*l./RIt,M 
501  CONTINUE 

WR I TF I  2  H  Ot ! I ),t*l,NI 
RFW1N0  ? 

RETURN 

FNO 


IIBFTC  SYMll 


SUBROUTINE  SYMIN  (S,N,IS1 
DIMENSION  S( 80, 1 ) 

N 1  *  N-l 
00  04  I  *  2,N 

II  *  I- 1  _ _ 

DO  04  J  *  1,11 
04  SI  I, J  1  =  0.0 

C  *  1 .O/SQRTI S (1,1)1 
S (  1 , 1  )  =  1.0 
DO  13  J* 1 ,  N 

n  sii,j»  =  sii,ji  *  c 

DO  21  K®2  ,N 

DO  17  J3  1  ,N 
K 1  =  K-l 
DO  17  1  =  1, K1 

17  SI K , J  1  =  SIK, J  )  -  SIIfKl  *  St  I • J ) 

IF(-S(K,K)  160,61,61 
60  CONTINUF 

C  *  1. O/SQRTI  SIK, KU 


no  m 

1*1, K1 

191 

SI  I,K  1 

=  0.0 

SIK ,K  1 

=  1.0 

DO  21 

J*  l  ,N 

21 

S I K  ,  J  1 

=  SIK.Jf 

DO  30  J*2,N 

J1  *  J-l 
.  DO  30  1=1, J1 

DO  30  K=J,fJ 

30  S(I,J)  *  S(I,J)  ♦  SIK, II  *  SIK, 01 
DO  35  J=  1  ,N1 
SIJ,J)  *  SIJ,J>  ** 2 
J  2  =  J*1 
DO  35  I*J  2»N 

3*;  S(J,J)  =  S(J,JI  «■  SI  I ,  J  J*+2 
S  , f N  1  *  S(N»N  1**2 
DO  4 7  1=1 ,N1 

12  «  U1 
DO  42  J*  1 2»N 

4?  S(J,I»  *  SII,J) 

GO  TO  62 

61  CONTINUE 
I S- 1 

62  CONTINUF 
RETURN 

END 


(ORIGIN 
(IBFTC  JACS2 


ALPHA 


SUBROUTINE  JACS3 

OINENSION  RI  160j80J  ,0(801 _  _  _ 

COMMON  P»NL»NtNF»LfKL»KKl,KK2l , NA ,El f EE , KK3L  f HH, KKK 
♦,NCfFFi,FF2,LLl 

*^Pn,QD,PW  _ _ _ 

*♦  NF1 
♦  ,  T IM 

♦.LI  _ 

♦,LLA#NP,LIB.LI6,LB,LE,JI 
00  53  1=1, N 

READ  (  2  I  ( Rl  I»  J  I , J*1 »N)  _  _ _ _ 

53  CONTINUE 
REWIND  ? 

06  N I  =  N*1 
061  N11=N-1 

07  N2  =  N*2 
08  00  10  I  =  N  If  N2 

09  00  10  J  =  l.N 

10  R(I,JJ  =  0. 

11  00  12  1  =  1,  N 

111  NI  =  N  ♦  I 

12  R(NI,  II  =  1  . 

R  I M  =  1  . 

DO  76  I  *  1, Nil 
II  =  IM 

13  00  75  L  =  1 , NL 
4B=0. 

DO  284  J=I1,N 
R  I J  =ABS  ( R  (  If  Jl  I 
AB=AMAX1( AH.RIJ) 
lP(P-RIJ140f42,42 
40  CONTINUE 
LR=  1 

0R=R< If  I  l-P(J,JI 
DRR=DR**2 

AK=SURT( DRR/I DRRf4,*R (I , J I ♦♦2 ) ) 

SO*  SI  ONI  l  •  ,0R  1 
A*SQRT( ( 1 .♦S04AK1/2.1 
2?  R  «  SORT  (  l »-A**2 ) 

221  C  =  SIGN  (  l.,R( I,J) I 
AC= A*C 
BC=B*C 

27  DO  252  K  =  i,N2 

U  *  RIK,  I  1*AC  ♦  R (K, J )*B 
R  I  K  ,  J  )  =  -R(K,ll*BC  ♦  R  I K  ,  J  I  *A 
252  RIK,  I)  =  U 

Ri  !,n*RU,n*AC+Rij,n*B 
R 1 J , J ) *-R I  I , J 1*BC *R  ( J  ,  J )*A 
RI 

RIJ, I  1*0. 

DO  283  K*1 , N 
R(IfKl«RIK,n 
R(J,K»*P(K,J1 
2R7  CONTINUE 
42  CONTINUE 


Ik.lk 


284  CONTINUE 

IFIP-AB  144,43,43 

44  CONTINUE 
35  CONTINUE 

43  CONTINUE 

IEIRIM-RI  I,H  145,46,46 

45  CONTINUE 
L  I*  l 

36  CONTINUE 

46  CONTINUE 

00  332  1*1, LI 
332  0(  1  )  =  R(t,n 

RETURN 
end 


^ORIGIN 
SIBFTC  JACSi 


ALPHA 


SUBROUTINE  JACS 

DIMENSION  R(160t80l ,D(80),A(80I 
* , DEI  80  I , OA ( 80 ) 

COMMON  P,NL,N,NF,L ,KL,KKl,KK2l , NA , E t , F E,KK3l , HH, KKK 
*,NC,FF1,FF?,LLI 
*,pd,qd,pw 

*,NF1 
*»  T IM 
*,L  I 

♦,LLA, NP,t IB,LIFtLBtLFf JI 
00  504  1  =  1, N 
READf  2MR(  I»J  )f  J*1«N) 

504  CONTINUE 
RFAD12MQEC  II  ,1  =  1,N1 
REWIND  ? 

PK=0. 

DO  61  I  =  1 f  N 
DA ( 11=1  .-DEI  I  1 

A  (  I  )=  l./SQRT(Pl)*DF  (  I|*QO*nA(!  »  ) 

I,K=AMAX1(PK»PW*0£(  I  >*A  (  I  )**2  I 
61  CONTINUE 
OO  6  3  T  =  1 f  N 
On  6?  J  =  1 »  N 

P(  I ,  J  I=AI  I  )*R1  I,J»*AU» 

6?  CONTINUE 

RIM  I  =  11  .-PW*DE(  I  II*AII  )**2*PK 
63  CONTINUF 

CALL  JACSIM(Rfn,P,N,NL,LII 
L=L  I 

DO  807  J=  1 1 L 
Of J  )  =  SQPTf  0 <  J l-PK  I 
DO  505  I*  1  ,N 
I N= T+N 

R( If J  )=RI  IN, J  *  *01 J  I 

505  CONTINUE 

WRITE  t4MRU,JJ,I  =  l,N  \ 

607  CONTINUE 
REWIND  A 
RETURN 
END 


xmm*  ^«uir  »;‘v 


il 8FTC  JACSIl 


SUBROUTINE  JACSIM  f R,D, P , N,NL »LI) 
DIMENSION  R< 160,801,01 801 
06  Ni  =  N+l 
061  Nll*N-l 
07  N2  *  N*2 
08  00  10  I=N  1,  N2 

09  00  10  J  =  1,N 

10  RU,J)  =  0. 

11  00  12  1=1, N 

111  NI  *  N  ♦  I 

12  RINI, I )  *  1. 

DO  36  I  =  l,Nll 
II  =  I+l 

13  00  35  L  =  l,NL 
AB=0. 

DO  284  J=  I  1 , N 
RIJ=ABS(R(  I,J)I 
AB=AMAX1I AB,RIJ) 

IE|P— RIJ)40,42,42 
40  CONTINUE 
IR=1 

OR=R  (  I ,  I  )-R  I J  ,  J  | 

DRR=DR**2 

AK=SQRT(DRR/IQRR44.*R  II,  J  1**2  M 
SO=SIGN( 1..0R) 

A=SQRT( (1 ,4SD*AKl/2.» 

22  8  =  SORT  I  1 A**2 I 
221  C  =  SIGN  (  l.,R(  I, J)} 

AC=A*C 

0C=B*C 

23  DO  252  K  =  1,\I2 

U  =  RIK,  I  ) *AC  ♦  R (K, J|*B 
RIK.JJ  =  -R(K,n*BC  ♦  RIK,JI*A 
252  RIK,  M  =  U 

r ( i , i  i=r i  i,n*AC4Pi j,  n*p 
R(J,JI=-Rf I,J»*BC4R(J,JJ*A 
R( I,J 1=0. 

RIJ, I )=0. 

00  283  K=1,N 
R I T,K)=R(K,I» 

RI J,K)=P(K,JI 

283  CONTINUE 

42  CONTINUF 

284  CONTINUE 
IFIP-AB  144,43,^3 

44  CONTINUE 
35  CONTINUE 

43  CONTINUE 

IF! L 1“ I  145,46,45 

45  CONTINUE 
36  CONTINUE 

46  CONTINUE 

00  332  1=1, LI 
332  01  II  «  PI  I, II 

RETURN 
END 


SQRIGIN 
SIBFTC  RARE  1 


ALPHA 


SUBROUTINE  RARE 

DIMENSION  Rl 80,80), AMI80,30), UMI80,30) ,WM!ll0,30) 

*t  D(  80 ) ,  U(  8 C  I 
*, A( 150) 

*,AA(150)  _  _ 

* »  DEI  80 ) , DA ( 80 ) 

*»  UEI 80 ) 

COMMON  P,NL,N,NF, L, KL  ,  KKL , KK2L , NA , E l , EF , KK3L ,HH, KKK 
*,NC,FFl,FF2,LLl 
* , PD» 00, PW 
* ,  NF 1 
*,  TIM 
*»  L  I 

*,LLA,NP,L  IB,LIE,LB,LE, JI 
T IM 1=  T  IM'*(  2 ) 

DO  701  I  ~  1  * N 
R  FAD( 2  )  ( R (  l,J ),J*l,N> 

70\  CONTINUE 
REWIND  2 
DO  702  J=1,L 
R  EADI 4 ) ( AMI  I , J )  ,  I  - l ,N  ) 

702  CONTINUE 
REWIND  4 
DO  63  1=1, N 
DEI  11  =  1. 

DO  61  J  =  1 , L 

DEI  I  )  =  !)E(  I  )*AM(  I ,  J  1**2 

61  CONTINUE 

D(  I )  =  1 ./SQKTI DE I  I  )  ) 

DO  62  J  =  1 , L 

AMI  I,J)=D(  I)*AM!I,J) 

62  CONTINUE 

dai  n=iDEi  ii-i.i/deiii 

DEI  1 )*1 . — C  A I  I  ) 

D|  I )=PO*DE| 1 1 +QO*DA( I  I 

63  CONTINUE 
L  N  =  L  *N 

DO  347  KKK *1 , KKl 

CALL  RARFDIR,  JM , A M,  WM ,D , N ,L  ,LN , U, KKK, KKL, AL 
*,EF 
*,UF 
*,C2 
*,  ALM 

*,PO,PW,QD,nA,OE) 

AKK-KKK 
A  A I KKK  )*ALM 
A I KKK  )*C2 

IFIP-AL ) 34 7,347, 3471 
347  CONTINUE 
3471  CONTINUF 

WR  !  TE  16,008  HAA  I  I  ),  1*1,  KKK) 

WRITE! 6,992  I 
<>92  FORMAT!//) 

WRITE(6,90eHA<  I),  1*1, KKK) 

908  FORMAT! 10F8. 4) 
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TIM 2* TIME!?) 

T  IM=T IM2-T IM1 
WRITE(8)C2,AKK,TIM 
WRITE(8HWWtI,ll,I  =  liLI> 

SI=0. 

DO  30  1=1, N 
D(  1 1  =  S1GN(  1., AMU  ,111 
SI=SI+Dl I  1 
30  CONTINUE 

SI=SIGN( 1. ,SI ) 

DO  34  1=1, N 
D( I|*0( I )*SI 
no  32  J= 1 , L 
AMI I,JI=AM(I,J»*S! 

32  CONTINUE 

WRITEIPIDI  n,OA(I»,DEm,fAM(I, J»#J=1,LU 
34  CONTINUE 

DO  1111  J* 1 ,L 
WRITE!  3KAMU,J1, 1*1, N1 
1111  CONTINUE 
2EWIND  3 
RETURN 
FND 


14-19 


UfiFTC  RARfPl 

SUBROUTINE  R AREDf P , UK , AM, WM,0 , N,L , IN.U t KKK ,KKL » AL 
*,F  P 
*,UF 
*,C2 
*,  ALM 

*,PO,PW,QO,CA,OE) 

0 {MENS  ION  R( 80,11 ,UMf  80, II ,AM( 80. 1 1 ,W*ftl(M  I ,D(80) ,UI80I 
* , OE ( 80  1 , DA ( 80  > 

*,UE<  80) 

00  815  [*t,N 
H  fit  1  )  =  1 .-PW*0E1 I 1 
nn  315  j  =  l »  l 

315  1)81  I,J)  =  I  AMI  !»J)  /Dill) 

DO  321  T=l»N 
IL  =  l  ♦  L 
on  .3?  I  J=  1  »L 
WMf  11  ,  J  1  =  0. 
on  321  K=  l»N 

■>21  WM(IL,J)  =  WMUI.J)  ♦  R(  I  »K  1  *  UMfK.Jl 
nn  3?8  1=1,1 
DO  328  J=I  ,L 

WM(  l , J  )  =  0.0 

on  327  K=1,N 

K.L  =  K  ♦  L 

32  7  W  M  f  T , J  >  -  WMtl.J)  ♦  WMfKL.ll  *  U*MK,J» 

328  W^IJ,!)-)  'I  l.J) 

r.i*o. 

on  ?  im.l 

C1=C1 *WM<  l  ,11 

?  continue 

00  336  K=1,L 

S  =  1  .  >/ SOH  T«W4(K,K>> 

(in  331  !*K,LN 

33\  WMU.Kl  *  WM1,K1  *  S 
K  l  =  K  ♦  1 

IF  (l-Kl  743,  343,  334 

3  34  OH  7  76  J  =K  1 ,  L 
on  736  !=J,IN 

536  W*MI.JI  -  v.Mfl.JI  -  WM  (  l ,  K I  ♦  WM(J,K» 

743  CONTINUE 

00  706  m,N 

net  ii*i. 

IL  *  l  *L 

00  701  J*1  ,L 

06(  I  1  *0F (  l  1-W81 It.J)**? 

701  CONTINOF 

JHFF-DH  I  II  fC4, 707,702 

702  CONTINUr 

00*  S18  T ( i l.-FFl/l 1.-PFI1III 

no  703  J*1 ,L 

W*t  IL , J  1  =WM( II , J 1 *00 

703  CONTI NUF 
DM  1»*M 

704  CONTINUE 

OA ( tl  =  l.-DF(  I  I 


UE(  II  =  OE(  I) 

D(  I  )  =  PD*DEimQD*DA<n 

ui i )*m 1 1 

705  CONTINUE 
Al*C2 
HH*0  • 

H1*U  l.-PW*0E  INI  I  /0IN»  )**2 
N1=N-1 

DO  11  I*l»Nl 

H1=H1*( < l.-PW*OEI 1 1 l/DI 1 1 1**2 
11*1*1 

DO  11  J*U,N 
H*R( T  *  J  ) 

11  HH*HH+H**2/ID( I»*D( J ) I 
HH*H1*2.*HH 
C2*Cl/HH 
ALM*0. 

DO  711  I  *  1  »N 
IL= I *L 

f)0  1 10  J*  l  »L 

ALM*AMAX1( AIM, ABS(AM( I , J 1-W" (I L , J )l » 
AMI  I «  J 1*WM I IL  t  J) 

710  CONTINUE 

711  CONTINUE 
AL*AHSI AL-C2I 
RETURN 

END 


SORIGIN 
SIBFTC  S IMP  1 


ALPHA 


SUBROUTINE  SIMP 

DIMENSION  A(80,30)*A1(80,30I,E(30,30) 
*,B(80,30l,S(30,30l,Ht30,30i,Dt80l 
*,01(80) 

*,02(80) 

*,00(80) 

*,  RC( 80, 30) 

*,G(80, 30) 

* ,  DL  (  30  ) 

*,CB( 30) 

*,  DF ( 30  ) 

COMMON  P,NL,N,NF,L,KL» KKL  »KK2L , NA,E1,EE,KK3L ,HH,KKK 
**NC»FF1,FF2,LLI 
*,PO,yD,pw 
*,NT  1 
* ,  T  IM 
*» L  I 

*,LIA,NP,LIB,LIE,IB,LE,JI 
M*L 
FMM=M 
FN  *N 
LIE*? 

LLE=4 
UE*3 
LLf *NL 
LL8«1 
U.F*0 
LL  1 
LL  F  =  1  0 
LLE  «  5 
LLc*?0 
ML  «  4 
ML«3 
Ml  *  1 
ML  *2 
NN«0 

DO  41  LLL*  l  »NF 1 
0')?  J*  1  , M 

»EAO(  3)(A(  I,J),t>),N) 

?.  CONTINUE 
REMIND  3 
DO  S’  I  *  l  »N 
00( l 1*0. 
no  51  J«1,M 
oo(  n«orm*A(i,j)**? 

51  CONTINUE 
DOm»SORTIOO(  III 
00  52  J*  1 » M 

I F INN 1 811,510, 511 

510  CONTINUE 

At  I , J  I  >  A ( I,J)/00l I  I 

511  CONTINUE 

9i i , j i*ai  i,:? 

52  CONTINUE 

53  CONTINUE 


Fl*0. 

00  42  LLL4=LIR,LLF 
ALL*Ltt4 
FN*llL4*ML 
F*2.*FM/t  2.*FM- l. ) 

FPI«F>1. 

FFF«l./|F-3.| 

DO  8?  J*l,M 
0F( J  |*0( J I 
00  82  1*1, N 
Gl  I,J  I  *BI  I  •  J  ) 

4?  CONTINUF 

»=1*0. 

00  20  Ll*l,Nt 
AU*Ll 

00  4  J*1,M 
Oil J»=0. 

01  J1*0. 

on  o  i  *  i ,  n 

OH  J  I  *01 1  J  )  *41  l ,  J  1**4 
01  J  1*01  J  | ♦  AOSI  dt  !,J))**FPl 
0  CONTINUF 

DIJ|*(01U)/0t  J|)**FFF 

01UI*0(JI 

00  4  I *J , N 

B(  I , J I *h(  1 1 J 1*01 J  > 

*  CONTINUE 
00  7  I*  1  ,N 
01 ! 1*0. 
no  A  J«l,N 
01  I  1*01 II+BI I , J|4*? 

6  CONTINUF 
on  ?  J»l,M 
All  t.JI*OI  IHAII.JI 
*CI It J • * AB S ( B (  ! , J I  |**F 
»  CONTINUE 

on  <»  i*i, n 

on  o  J*l,« 

M  l,J  1*0, 

Ft 

on  8  K  *  1 ,  N 

St  I,J  I* St  I  ,J|MIK,|  t*Al|K,J) 
Ft  If  J  l*rt  1  ,JI*AK,1  I'BCtK, J) 
8  CONTINUF 
4  CONTINUE 
CALI  SVPT 
on  10  !•!,* 

00  10  J*l  ,N 
H|  t ,  J  1*0, 

00  10  K*l,“ 

Ml  t , J  I *M|  I, JUSt  t,KI*rtK«JI 
10  CONTINUE 

c*o. 

Cl*0, 

no  240  J»l,M 

nt j i*o. 


*'•,***■•«• 


14-23 


on  23  !  =  1,N 
ou  i*oui*h<  i ,  j)**2 

21  CONTINUE 

D(J»*S0RT(0JJ»» 

C*C«-0U) 

00  24  1=1, M 
Hi  I  ,  J  l  =  H(  I 

24  CONTINUF 
3IJ|s0lJl/ClUl 
C1*CU0(J  I 

240  CONTINUF 

C=C/FMm 
C1=Cl/fMM 
no  26  J=l,M 
00  ">5  I  =  1 ,  N 
3t 1 ,J)=0. 

00  251  K*l,M 

Hi  1  •  J  1  *Bt  I  tJ»*M  !  »K|*HIK,J) 

?M  CONTINUE 

25  CONTINUE 
‘ >6  CONTINUE 

F  ?*F l 
F  1*C 

IFIP-I APSIEl-E 211 ♦?.» 52? ,528,528 
■»??  CONTINUE 
2  0  CONTINUE 
*>>4  CONTINUE 
LIA=UI  <■ 

ITf  mC,Cl*AL,Atl 
F2«FI 
*l  =  C 

IFIP-I ARSIEI-E21 l*2.ITC0»?08»708 
TOO  CONTINUE 

00  404  J»l»M 

08Ut*0. 

OO  403  I» l »N 

IF(  (M  l ,  J  I  ♦PI401 ,40?»40? 

401  CONTtNjF 

U4IJ)«CHI  JIM, 

40?  I.ONllNUf 
403  CONTINUE 

4*«AN|N|t A«,0*I JI  » 

40.  CONTINUF 

If  |  FM**-AW  l  264#  ?6t>»  40 
?*>*>  CONTINUT 
42  CONTINUE 
*0  CONTINUE 

IF!  AU-1*  I  702,  T06»?0? 

TO?  CONTINUE 

00  T04  JM,H 
01  J  )«0F  I  J  I 
00  T04  I  *  l  »N 
n(!,ji-gi i • j » 

?  14  CONTINUE 
704  CnNTINU* 


708  CONTINUE 

iF<NN)321»320t32i 

320  CONTINUE 

00  34  f  *i,r 
DO  34  J  =1,M 

8(i,ji=do!  n*au#ji 

34  CONTINUE 

321  CONTINUE 

00  36  I  *l,N 
WR  l  TE(  8  M  B  (  1 1 J  I » J=1  »LI  ) 
36  CONTINUE 

MRITEIBMOlJlf  J=l,LI» 
NN=  l 

41  CONTINUE 
RETURN 
END 


1U-25 


MBETC  SYMI2 

SUBROUTINE  SYKI3  IS»N) 

DIMENSION  SI  30, 30) 

38  N1  *  N-l 

01  DO  04  I  =  ?,N 

0?  II  =  1-1 

03  DO  04  J  *  1,11 

04  siiji  =  o. 

10  C  =  K/SQRTISI  1,1 1> 

11  SI  1,1)  =  1. 

12  DO  13  J=1,N 

13  Sli,J)  =  SU,J)  *  C 

14  DO  21  K=2» N 

15  DO  17  J=1,N 

151  K 1  =  K— 1 

1(5  DO  17  1=1,  K! 

17  S(K,J)  =  SIK,J|  -  SI  I » K )  *  $1  I , J) 
13  C  =  1  ,/SQRTI SIK.KI ) 

19  DO  191  I  =  1  ,K 1 
191  S ( I i K )  =  0. 

19?  SlK.K)  =  1. 

20  DO  21  J  =  1,N 

21  SI Kf J )  =  SIK.J)  *  C 

25  m  30  J* 2 ,N 

26  J1  =  J-l 

27  DO  30  1  =  1, J1 

29  no  30  K= J  , N 

30  S I  I ,  J  )  =  SII,J)  *■  SlK,n  *  S I  K ,  J  ) 

31  DO  35  J=  1 ,N1 

32  SIJ,J)  =  SI J, j)  **? 

33  J2  =  J*l 

34  DO  35  I=J2,N 

35  SIJ,J)  =  S|J,J)  f  S I ! , J ) *  *  2 
351  SIN.N)  *  S I N»  N ) **? 

3n  00  42  1= 1 ,N1 

40  12  =  I  ♦  l 

41  00  42  J  =  I  2  » N 

4?  SI J,  II  =  SIT, J) 

RETURN 

END 


14-26 


SDR  I GIN  ALPHA 

ilBFTC  OUPt II 

SU8ROUT INF  OUPL I 

DIMFNS  ION  A (80, 301 #D( 801, DA (80 1 ,DE(80) ,JI(80» 

*,01(80) 

COMMON  P, NL,N,NF,L  ,KL , KKL.KK2L , NA , E l ,E E ,KK3l ,HH, KKK 
*,NC,FF1 ,FF2,LL 1 
*,PD,QD,PW 
*,NF1 

*,tim 

* ,  L  I 

*,LLA,NP,L !8,L IE,Lfl,LE, JI 
NS  =  6 
NS=4 
NS  - 1 
NS=? 

HO  26  LS=l »NS 
L7*0 

00  24  LLL=l,NP 
WRITE (6, 999) 
cqq  CORMAU  lHt  ) 

00  22  LL I -L  IP, LIE 
00  20  LL=LR,l  E 
WRITE (6, 996) 

996  FORMAT (///) 

L  7=L  7  +  1 

REAC(8)LLL,N,LI,LLI,Ll,NFl 
o  F A0( 3 ) C2, AKK , T  IM 
kEAC(8M0(  I),  1*1, LI) 

WRITC(6,90I  )LLL,N,LI,LLI ,LL,NFl,C2, AKK  »TIM»(D(I)»I=l»LI) 

901  FORMAT ( 614, 3X , 3F7 .3, 3X , 7F 10. 3 ) 

WR l TE( 6, 997  ) 

997  FORMATMH  ) 

DO  2  1*1,  N 

READ!  8  1DH  I  )  ,  fOA  I  I  >  ,0E(  I  )  ,  (Ad  ,  J),J  =  1  ,L!  » 

2  CONTINUE 

I F( L !- 1 )8*  1 5,  8 
fl  CONTINUE 

LL  A  =  J I ( L  7 ) 

00  4  1*1, LLA 
READ(  8  IC,C  1 , AL  •  ALL 
WRITF(6,902  )C*C1,AL,ALL 

902  FORMAT ( 5F7 • 3) 

4  CONTINUF 

WR ITE ( 6 , 997 ) 

00  14  1*1,  N 
RFA0(8  )(0(J  I,J*1,L1  ) 

WRITE(6,905)(D(J)»J*1,LI) 

905  FORMAT! 78X,7F7.3) 

WRITE (6,906 )I,ni(I),OA(I),OF(I),(A(I,J),J*l,LI) 

906  FORMAT! 1H+, 1 3 , 1 X , 1 F4. 0, IX, 2F 7, 3 ,2 X , 7F 7. 3 ) 

14  CONTINUE 

WRITEI6.997) 

RFAD(8Kn(J)»J-l,LI  ) 

WRITE! 6,904 )(9(J|,J*l,LI) 

904  FORMAT! 2X, 7F7,’» 

GO  TO  18 


15  CONTINUE 

00  16  I*  1 » N 

WR!TE!6t9<nn,Dimf0AUJjDEIIJ,|AU_iJ)tJ»lili 
903  FORMAT!  13,  IX,  IF4.0,  1X,2F7.3,2X,  14F7.3I 

16  CONTINUE 

18  CONTINUE  .  _ 

20  CONTINUE 
22  CONTINUE 
24  CONTINUE 

WRITE(6,999I 
REWIND  8 

26  CONTINUE  .  _ _ 

RETURN 

END 
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